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Abstract 

We consider here the coexistence of first- and third-order integrals of 
motion in two dimensional classical and quantum mechanics. We find 
explicitly all potentials that admit such integrals, and all their integrals. 
Quantum superintegrable systems are found that have no classical analog, 
i.e. the potentials are proportional to h , so their classical limit is free 
motion. 



1 Introduction 



In classical mechanics, an n-dimensional Hamiltonian system is called Liouville 
integrable if it allows n functionally independent integrals of motion in involu- 
tion (including the Hamiltonian), that is 



{H,.Y,) = 0. 
{X,,X j } = 0,Vi,j. 

The Hamiltonian H = H (xi, ...x n ,p\, ...,p„) and the integrals of motion 
Xi = Xi(xi, ...x n ,pi, ...,p n ) must be well defined functions on phase space (|^|, 
|l4j]). The system is superintegrable if it allows more than n functionally in- 
dependent integrals, n of them in involution. The best known superintegrable 
systems in n dimensions are the harmonic oscillator V = lot 2 and the Coulomb 
potential V — f, both of them allowing 2n — l independent integrals of motion, 
the maximal number possible for an interacting system. Bertrand's theorem 
((|, D) tells us that these are the only rotationally invariant systems for which 
all finite trajectories are closed, a fact intimately related to their maximal su- 
perintegrability. 

In quantum mechanics, a Hamiltonian system is said to be integrable if there 
exists a set {Xi} of n well defined, algebraically independent operators (includ- 
ing the Hamiltonian) that commute pairwise. It is superintegrable if it possesses 
further independent operators, {Yj} that commute with the Hamiltonian. The 
Yj do not necessarily commute with each other, nor with the X^. 

The definition of the independence of quantum operators is not unique, and 
this may give rise to different types of quantum superintegrability. 

A good working definition, which may be appropriate for applications in 
quantum mechanics, soliton theory and for instance in the study of the Huygens 
principle, is that operators are considered independent unless one of them can be 
expressed as a polynomial in the others ( |[ [f6|, ) . The fact that commuting 
operators can be useful even if they are functionally dependent in the classical 
limit was clearly demonstrated by Hietarinta (jl6j [l?], [l8|, This definition 

is in itself not quite satisfactory since it ignores more general polynomial or 
functional relations between integrals. This may lead to important differences 
between classical and quantum integrability. Moreover, it is not appropriate 
for nonpolynomial integrals. Finding an appropriate and rigorous definition 
of the independence of quantum operators is not an easy problem, but it is 
worth investigating as wrong or ambiguous definitions may give rise to incorrect 
results. For a discussion of related problems, see e.g. U) and @ 

Previous systematic searches for superintegrable systems concentrated on 
integrals of motion of at most second order in momenta (|| ||, |l0[ |l3|, ^J, |3l||). 
This "quadratic superintegrability" has been shown to be related to multiscp- 
arability of the Schroedinger or Hamilton- Jacobi equations. More recently, it 
was related to generalized symmetries (^8|) and exact solvability ([p9|). 

Quadratic superintegrability has been considered in spaces of nonzero con- 
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stant curvature (p3l |2(|) and of nonconstant curvature (|2~]]]) . For superinte- 
grable systems in n dimensions see (|27[). 

The purpose of this article is to start a systematic search for superintegrable 
systems with higher order integrals of motion. We consider a two-dimensional 
real Euclidian space with a one-particle Hamiltonian; 

H=l(pl+P 2 y )+V(x,y). 

We request the existence of two additional integrals of motion, one of first 
order in the momenta and the other of third order. 

The classical and quantum mechanical cases will be treated separately. When 
second order integrals of motion are considered, classical and quantum inte- 
grable and superintegrable potentials coincide. For third order integrals this is 
no longer the case (as was pointed out by Hietarinta in [jlfj)). For integrable 
systems with third or higher order integrals in classical mechanics, see also 

§ n n [§ II n 

2 Conditions for the existence of a third order 
invariant in classical mechanics 

We are looking for a classical integral of motion that is a polynomial in the 
momenta with coefficients depending on the spatial coordinates, i.e. 



that Poisson-commutes with the Hamiltonian; 

= {H,X}, 

„ Pi +Pl _l rr, \ 

H = — 2 — +V{x,y). 

We can simplify our search by using the fact that equation (| 
X is a constant over any trajectory: 

dX dX . dX . 

at aqi dp l 



with 



Pi = -V qi (qi,q 2 ), 
<l% =Pi- 



(2.1) 
implies that 

(2.2) 
(2.3) 



If we write explicitly X in (2.2), we find 
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E( dfjk j+i k 
{-to* P > 

j+k=X V 



dfjk j k+i 



fjkV x jp[ 1 pl - fjkV v kp\pl 



0. 

(2.4) 



Since the monomials pipe's form a basis, the coefficients for each (a, b) must 
vanish separately, thus (2.4) gives relations between the with odd and even 
i + j separately. If we are looking for an integral of odd (even) degree in the 
momenta, the even (odd) terms will play no role and we can without loss of 
generality consider only integrals that have terms only of odd (even) parity. 
Moreover, we may notice, in (2.4), that the terms of leading order in the p^s 
imply a relation independent of V between the fij with i + j = n. This allows 
us to find immediately the form of the leading order terms, so the integral of 
motion in the third-order case takes the form 



X= ^2 A ljk p\p : ' 2 L k +g 1 (x,y)pi+g 2 (x,y)p 2 , 

i+j+k=3 

L = xp 2 - ypi. 
where the are arbitrary real constants. 



(2.5) 



The requirement =£■ — and the Hamilton equations ((2.3)) yield four 
equations 





(gi)x 

(S2)y 
(gi)y + (ff2)x 



9iV x +g 2 V y , 
3fi(y)V x + h(x,y)V y , 
h(x,y)V x + 3U(x)V y , 
2{h{x,y)V x + h{x iy )V y ) 



(2.6) 
(2.7) 
(2.8) 
(2.9) 



where 



hiv) 

h(%,y) 
h{x) 



-^3ooy 3 + Mwy 2 - Ai 20 y + A 030 , 
3A 300 xy 2 - 2A 2W xy + A 2Q1 y 2 + A 120 x - Amy + A 2i, 
-3A 30Q x 2 y + A 2W x 2 - 2A 201 xy + A nl x - A W2 y + A ou 
A 30 qx 3 + A 20 ix 2 + A W2 x + A 003 , ■ 



Requiring that equations (2.7),( |2.8| ), and ( |2.9| ) be compatible, we obtain a 
linear compatibility condition for the potential, namely 



— — f?,V XXX + (2/2 — 3/4) V X Xy + ( — 3/l + 2/3) V X yy — f 2 Vyyy 

+ 2 (f 2y — fax) V xx + 2 (—3 fiy + f 2x + f-3 V — 3/4x) V xy + 2 (—f 2y + f 3x ) V yy 



+ i~^flyy + 2/2x1, — f'ixx) V x + ( — f 2 yy + 2f 3xy — 3f4 XX ) Vy 



(2.10) 



4 



Requiring that all four equations ( |2.6| ) (2.7),p^),(^S|) be compatible, we 
obtain further third-order equations for the potential, this time nonlinear ones. 
They are limit case (for h — ► 0) of the corresponding quantum compatibility 



conditions (3.7) to (|3.9|) given below 



These conditions, together with (2.1C), form an overdetermined system for 



the potential V(x,y). The solution space will hence be rather restricted. In- 
deed, in 1935, Drach (JjJ) posed the problem of finding classical Hamiltonian 
systems with one third-order integral. In a complex Euclidian space i?2(C) 
he found 10 such potentials, each one depending on arbitrary constants, not 
however on arbitrary functions. We recall that in the case of second order in- 
tegrals, one obtains four families of potentials, each of them depending on two 
arbitrary functions of one variable [Si]). They are the four most general 
potentials that allow separation of variables in cartesian, polar, parabolic and 
elliptic coordinates, respectively. 

3 Conditions for the existence of a third order 
invariant in quantum mechanics 

Here we are interested in the existence of third-order operators,i.e. 

3 

x = J2 p ij( x >y)p\pi> 

i+j=Q 

Pi = ~ihd x , p 2 = -ihdy, 

that commute with the Hamiltonian. An equivalent way of writing this operator 
is 

3 

X= {P^(x,y), P \pi}. 
i+i=o 

Here the bracket means the anticommutator: 

Each of these anticommutators can be expressed as 

{f,p{p k 2 } + + i{f,pb%}- = r>-' . /' -/'J /<U + i{^m[f],M}- 
Hence we can write the operator X in the form 

X = X+ + iX~, 

where X + and X~ are self-adjoint operators. As the Hamiltonian itself is self- 
adjoint, X> = X + — iX~ must also commute, as well as X + and X~ . These 



5 



two last operators commute under the same conditions, so we may restrict our 
search without loss of generality to self-adjoint operators. This turns out to be 
quite useful in view of the following result 

Proposition 3.1. For each self-adjoint integral of motion of order n, there 
exists one integral of order n with definite parity, i.e. 

[f ] n-2j 

x - = E E {Pn^A^y),p k iP n 2~ 23 ~ k }, (3.i) 

j=0 k=Q 

where P is a real function. 

Proof. This is simply due to the fact that we have a real Hamiltonian and 
a purely imaginary momentum operator, so terms of even order, which are 
real, must commute independently of the terms of odd order, which are purely 
imaginary. 

□ 

In the case n = 3 we restrict ourselves to third-order integrals of the form 
^3 = E {fij( x >y)>PiP2} + {9i(x,y),Pi} + {92(x,y),P2}- 

Requesting 

0=[H,X], 

H=^(pl+Pt)+V(x,y), 

we find a set of 15 differential equations, of which the first nine can be explicitly 
solved to give 



P1P2} + {9i(x,y),Pi} + {g2(x,y), P 2}. (3.2) 

i,j,k 
i+j+k=3 

where the Aijk are arbitrary real constants. So far this is similar to the classical 
case. 

Remark 1. The argument used in demonstrating proposition \S. i| can be gener- 
alized to any expression involving the anticommutators of self-adjoint operators 
homogeneous in the pt 's, for example to terms of the form {L\ : p\p^} , as long 
as the coefficients of the pi 's are real. 
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We could get rid of the h and m factors by a dilation of the undetermined 
functions, 

h 2 - 

V(x,y) = 2^ V ( X 'V)> 
gi{x,y) = h 2 g[(x,y), 
92{x,y) = h 2 g' 2 {x,y). 

This is equivalent to setting fi and m equal to one, which we could do, but 
we prefer to keep track of the dependence on Ti (while setting m = 1), in order 
to see the classical limit. 

We are left with a set of 6 equations, two of which are consequences of the 
other four, so that, as in the classical case, we have to solve four equations: 



= g{V x + g 2 V y -(fiV xxx + f 2 V xxy + hV xyv + fiV yyy 

' ' (3.3) 



+ 8^300^^ - yV x ) + 2 (A 2W V X + A 20 iV y ) ) , 

(gi)x = 3fi(y)V x + f 2 (x, y)V v = hi, 
{g2) v = h{x,y)V x + 3f4(x)Vy = h 2 , 
(gi)v + (.92). - 2 (f 2 (x, y)V x + f 3 (x, y)V v ) = h 3 . 



(3.4) 
(3.5) 
(3.6) 



E qua tions (3.4) to (|3.6| ) are t he s ame as in the classical case, however equa- 
tion ( |3.3| ) differs from equation (2J3) by the terms proportional to Ti 2 . Both in 
the classical and quantum cases we can eliminate <?i and g 2 and obtain compat- 
ibility conditions for the potentials. 

We shall write these in a unified manner for both cases. One such compat- 
ibility condition is the third-order linear equation ( 2.1 0| ) . To write three more 
conditions we introduce the notation 




and use hi,h% and /13 introduced above. In the classical case we have <p2 = 0. 
The three (independent) nonlinear compatibility conditions are 

</>i {h 3 (j)i + h 2 tfil + (j>i(j) 2y + (j) 2x + hi) \ 
<Plx + (Pl<Ply J 
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/ (j)\h b + fafay + 01^6 + fax + hj 
V <t>lx + fafay 

hi (faxy + <fly) + ^5 (<\>\faxy ~ fa x ~ ^fafaxfay) 
+ K (fafaxy - faxfay) - (hiy + fah 5x ) (01a; + fafay) 

= -02x (01xj/ + 4>ly) + 4>2y {faxfay ~ faxy fa ) + faxy {fax + fafay) ■ 

(3.9) 

In the quantum case these are fifth order equations for the potential. They 
can be used to express faxy, faxx and fa yv in terms of 02a;, fay an d fa- In the 
classical case we set fa — 0, but the equations remain independent. They can 
be used to determine fa xy , faxx and fa yy in terms of fa x , fay and fa. The 
nonlinear equations for V(x,y) are third order ones in the classical case. 

In deriving these equations we have assumed 

fax + fafay + 0, 
faxy fa ~ fax fay 7^ 0. 

The cases when the above conditions do not hold must be considered sepa- 
rately. This will actually be the case for potentials considered in this article. 

We also mention the interesting fact, already noticed by Hietarinta (|l7|]),that 
a classical integrable potential is also quantum integrable, if and only if it re- 
spects the compatibility condition 




f\Vxxx + J'lVxxy + /3^a 



xyy 



+ hVyyy + 8A 30Q (xVy - yV x ) + 2 {A 2W V X + A 2m V y ) = 0. 

(3.11) 



In that case the equations are invariant under a simultaneous dilation of the 
potential and t he g^s. Thus any potential that is a solution to both (2.6) to 
(U) and ( |3~li"| can be multiplied by an arbitrary factor, which can be used to 



"absorb" the h factor so the solution does not vanish in the classical limit 

Even if a classical superintegrable potential does not satisfy this relation, 
there could exist corresponding quantum superintegrable systems. In that case, 
though, the equations are not invariant under a dilation of the potential as in 
the previous case, so terms that do not satisfy both (3.11) and (2.6) to (2.9) 
must be proportional to h 2 and vanish in the classical limit. 

We will show that condition ( 3.1l| ) cannot be the consequence of equations 
(|Q|), ( |3.5| ) and (3.6) for in that case the classical and quantum integrable po- 
tentials would be the same. 



4 Superintegrable systems with one third order 
and one first order integral 
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4.1 Integral of first-order 



A potential V(x, y) allows an integral that is of first order in the momenta if and 
only if it is invariant under either rotations or translations. Thus the potential 
must satisfy 

aL 3 V + bpxV + cp 2 V = 0. 

Without loss of generality, we can take the potential and first order integral 
to be one of the following: 

t a^O: V = V(r), X = L 3 

• a = 0, b 2 + c 2 ^ 0: V = V{x), X = p 2 . 



4.2 



Quantum and classical superintegrable potentials in- 
variant under rotations. 



Compatibility conditions obtained from equations (2.6) to (2.9) or (3.3) to (3.6) 
leave us with only two possibilities, namely 



V 



a 



V = co 2 r 2 



The Coulomb potential and the harmonic oscillator, which are the best- 
known superintegrable potentials in any dimension. In addition to angular 
momentum L3, the Coulomb potential in E 2 allows two second order integrals, 
namely the components of the Laplace-Runge-Lenz vector: 



X? = {L 3 , P1 }- 
XC = {L 3 , P2 } + 



2ay . 

r 
2ax 



The harmonic oscillator, in addition to angular momentum, allows two sec- 
ond order integrals which are the components of a quadrupole tensor; 



vh 1„2 1 2 , 22 22 
X[ = --pi + -p 2 +uj x -uj y ■ 

x 2 = -P1P2 + 2uj 2 xy. 

Commuting (or Poisson commuting) second order integrals, we in general 
find third order integrals. 

The third order integrals obtained for these potentials are indeed direct 
consequences of integrals at order one and two. 
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4.3 Classical superintegrable potentials invariant under 
translation. 

In the classical case, the remaining equations are readily solved. If we set V y = 0, 
equations (j2.6|) to Q2.9D simplify to 



= 5i 

= A 30 o = A 2W = A 120 = A 030 
(92) y = h{x,y)V x 
{92)x = 2f 2 (x,y)V x 

We can at once set Aq 2 i and Aoq 3 to 0, for they correspond to trivia l co n 
stants of motion, p 2 and Hp 2 , that can be subtracted from the constant (2.5). 



The compatibility condition between the two last equations forces one of the 
three following conditions to be satisfied (up to a translation in x). 

V = ax (4.1) 
V=- 2 (4.2) 



x 



1201 



Am = A W2 = A012 = (4.3) 



The first two potentials correspond to superintegrable systems that have one 
first and at least one second order integral. Their third order integrals can be 
obtained by commutation of these lower-order ones. 

The last conditions forbids the existence of a nontrivial third-order commut- 



ing operator for any other potentials than (4.1) and (4.2) 



4.4 Quantum superintegrable potentials invariant under 
translation. 



Here the situation is more interesting. Equations ( p.3[ ) to ( p.6| ) reduce to 

= giV x - — (fiV xxx - SyA 300 V x + 2A 210 V X ^ (4.4) 

(91) x = 3fi(y)V x (4.5) 

(92) y = h(x,y)V x (4.6) 
(9i)y + (92) x = 2(f 2 (x,y)V x ) (4.7) 

The linear compatibility condition leads to two equations (since coefficients 
of y° and of y 1 must vanish separately), namely 

= (A 210 x 2 + A in x + A 012 )V XXX + A{2A 210 x + Am)V xx + I2A210V;; 

= (3^ 300 a: 2 + 2A 201 x + A W2 )V XXX + i(6A 30Q x + 2A 201 )V XX + 36^300^. 

(4.8) 
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The two equations are similar and easy to solve, but it turns out their only 



solutions that also satisfy (4.4) to (4.7) are again the potentials V = ax and 
V = a/x 2 . Their third-order integrals in general are direct consequences of 
lower-order commuting operators, that is they can be obtained by commuting 
their second order integrals. In the V = -% case, we find three third order 
integrals, 

X 1 = {Llp 2 } + a{2^,p 2 } 

y 

X 2 = {L 3 ,PiP2} - a{4— ,p 2 } 

^3 =pjp2 - a{-^,P2} 

The integrals X 2 and X3 can be obtained by commuting Xi with the first-order 
integral p 2 . 

In the particular case V = jy, we find four more integrals, again related to 
each other by commutation with p 2 ; 



Ti 2 6v 2 h 2 — 3w 3 

Xi = L\ + ^r{— + 2x, P2 } + -Tri^f- ~ 

x 5 = {L 2 , Pl } - h 2 {^, P2 } + '4-{^L. + 



x& = {l 3 ,pI} - h 2 {-,p 2 } + h 2 {— ^,pi} 



x 7 =pl + Y {-, Pl } 

In this case we find nine third order integrals, two of which are trivial (Hp 2 
and p 2 ), an d four are purely quantum integrals. In the classical limit they 
correspond to integrals of the free motion. Only the first three can be associated 
with the corresponding classical integrals of V = -% . 

Th e most interesting potentials are obtained by setting all the Aijk involved 
m JO) equal to 0. The expressions for /i,/2,/3,/4 greatly simplify and equations 



(4.4) to (4/7) can be solved directly. The nonlinear compatibility condition for 



these four equations reduces to 



h 2 V'(x) 2 = W{xf + aV(x) 2 + (3V{x) + 7, (4.9) 



where the a,/3,7 are arbitrary real integration constants. Equation ( |4.9| ) is the 
well-known equation for elliptic functions which can be rewritten as 



h 2 V'(x) 2 = A{V{x) - Ai){V(x) - A 2 )(V(x) - A 3 ). (4.10) 
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The constants Ai are either all real, or one of them is real and the other two 
are complex conjugated. If all three constants are real, we obtain either finite 
or singular potentials of the form 

V\ = (hu) 2 k 2 sn 2 (Lux,k), 

= (M 2 ( 4 - n ) 

sn 2 {<jjx, k) ' 

respectively. 

If we have e.g. ^3 = A\ and ImA^ ^ 0, we obtain the singular potential 

y (^) 2 

3 2(cti(lux, k) + 1) 

(throughout we have < fc < 1, w 6 R). 

The special cases with k = or k = 1, which arise when two roots coincide, 
can be expressed in terms of elementary functions. The most interesting example 
is the "soliton" potential, 



Vla = (M 2 

cosh 2 (lux) ' 

obtained by setting k = 1 in V\. If we set k — 0, or k = 1 in V2, we get a 
singular periodic, or nonpcriodic potential, respectively, namely 



V2a- ^ 



(Km) 2 
sin (cox) 
(Hlu) 2 



sinn {lux) 

For all these potentials lu is an arbitrary constant, hence there exist potentials 
of arbitrary amplitude for all nonzero values of U. 

Finally, if all roots coincide, we reobtain the known supcrintegrable potential 



y t -K 



X 2 ' 



which explains the extra integrals found previously for that potential. 
The other potentials Vi,V2,V^ also satisfy 



— — W = a, a = A 1 +A 2 + A 3 , 

4 V x 
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a consequence of (4.9) 



The two nontrivial integrals of motion for all these potentials can be written 

as 



Xi = {L 3 ,P 2 i} + {(a - W(x))y, Pl } + {-ax + 2xV{x) + / V{x)dx,p 2 } 



X 2 =pl + ±{W(x)-a, Pl }. 



(4.12) 



The second integral can be trivially obtained by the commutation of the first 
one with p 2 - 



5 Conclusion 

We have found all potentials in two-dimensional Euclidian space E% that allow 
one first- and at least one third-order integral of motion. In the classical case 
the result provides no new superintegrable potentials; all the potentials found 
allow second order integrals and the third order integrals are consequences of 
the second order ones. In the case of quantum mechanics the result is quite 



different. Any potential satisfying the elliptic function equation (4.9) will be 
superintegrable in the above sense, i.e. will allow the first-order integral p 2 and 
two nontrivial third order integrals. All those "behave well" in the classical 
limit, that is they are proportional to fi 2 and therefore their classical limit is 
the (superintegrable) free motion. 

No new superintegrable systems are found for rotationally invariant poten- 
tials V(r), neither in the classical, nor in the quantum case. Thus all potentials 
found above are of the form V = V(x), i.e. are actually one-dimensional. The 
problem however remains two-dimensional as the kinetic energy and the inte- 
grals of motion also involve the y direction. 

There is also an interesting link with soliton theory (pi). All new super- 
integrable potentials obtained above are also translationally invariant solutions 
of the Korteweg-de Vries equation. The same potentials occur in the rational, 
trigonometric and elliptic Calogero-Moser-Sutherland models (||). 

The difference between classical and quantum integrable and superintegrable 
systems with higher order symmetries makes the systematic search for such sys- 
tems very interesting. First of all, Drach's study of classical integrable systems 
should be completed. His systems are really complex ones and most of them do 
not exist in real Euclidian space. Moreover it is not clear how complete his list 
is. On the other hand, Rahada (|25|]) has shown that 7 out of 10 Drach poten- 
tials are "reducible" in the sense that they are superintegrable and allow two 
second-order integrals. The third order integral found by Drach is the Poisson 
commutator of the second-order ones. 

The problem of classifying quantum systems with third-order integrals re- 
mains open and the conditions of Section |^ provide the means for finding all 
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such systems. 

Work is in progress on superintegrable systems in two-dimensional Euclidian 
space with one second and one third-order invariant, as well as with two third 
order ones. 
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